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Abstract. We present atomic-scale computer simulations in equiatomic Llo-CoPt where Molecular Dynam- 
ics and Monte Carlo techniques have both been applied to study the vacancy-atom exchange and kinetics 
relaxation. The atomic potential is determined using a Tight-Binding formalism within the Second-Moment 
Approximation. It is used to evaluate the different saddle-point energies involved in a vacancy-atom ex- 
change between nearest-neighbour sites. The potential and the saddle-point energies have been used to 
simulate the relaxation of the long-range order in CoPt using a Monte Carlo technique. A vacancy migration 
energy of 0.73 ± 0.15 eV and an order-disorder transition temperature of 935 K have been found. 



PACS. 61.43.Bn Structural modeling: computer simulation 
66.30.Fq Selfdiffusion in metals and alloys 



64.60.Cn Order-disorder transformations 



1 Introduction 

Ordering kinetics in intermetallic compounds have been 
the topic of many experimental and theoretical works. Sys- 
tems presenting a phase diagram derived from the Au-Cu 
canonical phase digram are among the most studied. Ow- 
ing to their high magnetic anisotropy, some of them in 
the iron-group metals and platinum group metals are of 
outstanding technological importance 1,2,3,4,5 . A good 
knowledge of the ordering process and its dynamics is thus 
a necessary step in any extensive research on these sys- 
tems. Among the energetic parameters that drives diffu- 
sion, the saddle-point energies and the migration energies 
are generally less known and hard to measure in ordered 
intermetallic compounds. In pure metals and random al- 
loys, the migration energy can be deduced, for example, 
from stage III of resistivity recovery during annealing af- 
' ter low-temperature irradiation |HJ or thorough analysis of 
residual resistometry along isothermal and isochronal an- 
nealing series Those two methods are, however, very 
sensitive to the microstructure of the samples and to any 
impurities or defects. Based on the earlier work of Flynn 
0, Schober et al. J0| have proposed a model for deter- 
mining the migration energy from the phonon density of 
states in pure fee and bec metals. This model has been 
extended latter on to AB3 compounds with LI2 struc- 
ture jllj . A prerequisite for such a determination is the 
knowledge of the phonon dispersion at different tempera- 
tures and states of order. In parallel to these experimen- 
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tal approaches, computer simulations, based on Molecular 
Dynamics (MD) and Monte Carlo (MC) techniques, allow 
one a wide range of kinetics relaxation problems and ther- 
modynamic properties to be studied. Offering a possibil- 
ity for investigating atomic migration in terms of crystal 
energetics treated by means of advanced solid-state theo- 
ries, MD simulations are, however, technically limited to 
rather small samples and short periods of real time. Con- 
sequently, being very powerful when considering the dy- 
namics of single atoms in a solid ^5|> the method is less 
useful for studying net kinetic effects such as, e.g., long- 
range order (LRO) relaxation. Most of simulation studies 
of structural kinetics are thus carried out by means of 
the MC technique (for references see numerous works of 
Binder, e.g., Ref. Crucial for the success of any sim- 

ulation, the interatomic potential can be deduced from 
first-principle MD, providing an accurate description of 
the atomic interactions, but requiring enormous computa- 
tional time and a limited number of particles. To a great 
extent, this limitation can be overcome by using empiri- 
cal or semi-empirical potentials which have the advantage 
of reproducing fast and with satisfying accuracy the ther- 
modynamic and structural properties of materials. Very 
satisfying results have been obtained, in this way, in tran- 
sition metals and alloys |14U15j . 

The Co-Pt system has interested many researchers from 
a practical viewpoint as providing catalyst materials |16| , 
but from a basic viewpoint as a magnetic system with 
a coupling between chemical and magnetic ordering [171 
I18H19| . Around the 50/50 stoichiometry, CoPt orders in 
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the Llo structure, made of alternating pure cobalt and 
platinum (100) planes |2()| . This very anisotropic chemi- 
cal order is accompanied by a strong magnetic anisotropy 
and a tetragonality (c/a < 1). The pronounced anisotropic 
properties of CoPt are at the core of the present renewed 
interest in this system |21II22II23| . In thin films, Llo-CoPt 
alloys display, in addition to a high magnetic anisotropy, 
a magnetisation that is perpendicular to the film surface 
, thus making such films good candidates for magneto- 
optical-storage devices [21]. This system has been studied 
extensively to understand the origin of the asymmetry in 
its phase diagram. Various explanations have been pro- 
posed: for example, an effect of many-body interactions 
|19|. an influence of magnetism, and variations of atomic 
interaction with composition J22122J- Calculation of this 
phase diagram may appear quite simple a priori, as it is 
based on the fee lattice at almost all temperatures and 
compositions, but it still remains a challenge today. The 
kinetics of atomic ordering in the ordered phases in this 
system has also been studied. The activation energy of 
changes in chemical LRO has been determined by resis- 
tivity measurements to be 3.18 eV and 2.52 eV for CoPt3 
|27] and CoPt gH] respectively. 

It is well known, that the vacancy-driven mechanism, 
namely, the vacancy-atom jump between nearest-neighbours 
sites, is the dominant microscopic process of diffusion in 
dense phases |29| . During such a jump, the energy bar- 
rier formed by the nearest-neighbours that the jumping 
atom has to overcome for completing the jump plays an 
important role. In this paper, we present MD simulations 
based on the Tight-Binding Second Moment Approxima- 
tion (TB-SMA) formalism |3U], which is well adapted to 
transition metals and alloys [HJ, to calculate the saddle- 
point energies for the different kinds of atomic jumps. The 
calculated saddle-point energies will be used, in conjunc- 
tion to the potentials deduced from the TB-SMA, to simu- 
late the LRO relaxation in CoPt and deduce the migration 
energy using a MC technique. 



2 Saddle-point energies calculation 

In this section, we present the many-body potential and 
its parametrisation procedure. Once well defined, the po- 
tential will be used to simulate the atom-vacancy jump 
mechanism and deduce the energy barriers for the differ- 
ent kinds of atomic jumps to a nearest-neighbour vacancy 
in the Llo-CoPt structure in the case of a rigid and a 
relaxed lattice. 

The potential we have used is based on the approach 
proposed by Rosato, Guillope, and Legrand [I3221GG3 where 
the energy Ei of an atom at site i, derived in the TB-SMA 
formalism |30j . is written as the sum of two terms; an at- 
tractive band energy (-£,•) and a repulsive pair interac- 
tion term (Ef). The band term is obtained by integrating 
the local density of states up to the Fermi level this 
gives rise to the many-body character of the potential nec- 
essary to account for surface relaxations and reconstruc- 
tions whereas the repulsive term is described by a sum of 



Born-Mayer ion-ion repulsions. When replacing the realis- 
tic density of states by a schematic rectangular one having 
the same second- moment |15j . one obtains: 




and 

ET= Ajjexpl-pui^-l)} (2) 

I and J indicate the chemical species Co or Pt. £jj is 
an effective hopping integral and q/j describes its de- 
pendence on the relative interaction distance, p/j is re- 
lated to the bulk modulus of the alloy under consider- 
ation. In our case r^ 1 is the first-neighbour distance in 
the metal I. Tq J as a free parameter and it is different 
from the first-neighbour distance rg =( r^-l- r o' 7 )/2- r ij 
is the distance between atoms at sites i and j. E\ and 
E\ are cancelled beyond a cut-off radius r c which has 
been chosen as the second-neighbour distance. From r c 
up to the (Co-Co) third-neighbour distance, the poten- 
tial is linked up to zero with a fifth-order polynomial in 
order to avoid discontinuities, both in the energy and in 
the forces. The TB-SMA potential applied to intermetal- 
lic compounds gives very reliable results with interactions 
up to the second nearest-neighbour as long as interme- 
diate and low-temperature ranges are considered. Long- 
range interactions are required in order to reproduce high- 
temperature properties which are usually beyond the ca- 
pability of short-range potentials. Considering the tem- 
perature range (600-980 K) and the goal of our simula- 
tions, the approximation of interactions up to the second 
nearest-neighbours is rather justified. Furthermore, x-ray 
and neutron diffuse scattering measurements have been 
used in conjunction with inverse cluster variation method 
to calculate the effective pair interactions in Ll2-Co3Pt 
and in Ll2-CoPt3 [251126] . In both systems the first and 
second interactions are the predominant ones. The same 
calculations are now in progress in equiatomic CoPt |35j . 

The potential parameters A/j, p/j and q/j are a 
priori unknown and will be determined by fitting the po- 
tential to the universal equation of state driving the vari- 
ation of the potential with distance |14| . This procedure 
requires usually the knowledge of the cohesive energy, lat- 
tice parameters and bulk modulus of the system. We have 
taken these parameters from the literature ■ We aimed 
during the fit to reproduce correctly the formation ener- 
gies and the lattice parameters of the ordered Llo and the 
disordered Ai phases of CoPt . 

The final parameters of the potential are summarised 
in Table I. In Table, II we report together the calculated 
lattice parameters and the formation energies and com- 
pare them to the experimental data. A good agreement is 
achieved for both ordered Llo and disordered Al phases. 
We have attempted to parameterise the TB-SMA poten- 
tial for the Ll2-CoPt3 phase but a discrepancy of more 
than 40% appeared between the experimental and the cal- 
culated formation energies. For this reason, further calcu- 
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lations, aiming at the determination of the saddle-point 
energies, have been restricted to the Llo phase. 

In order to calculate the saddle-point energies, we have 
used a simulation box containing 10 3 Ll fcc-based cells. 
In accord to the Llo structure, there are as many Co as 
Pt atoms arranged on (100) alternating crystallographic 
planes of pure Co and Pt using periodic boundary condi- 
tions in the three space directions. The total energy of the 
system, containing a single vacancy, was monitored dur- 
ing the four kinds of atomic jumps to a nearest-neighbour 
vacancy in the case of a rigid and a relaxed lattice: Co 
to Co-site vacancy (Co— >Co(V)), Co to Pt-site vacancy 
(Co->Pt(V)), Pt to Pt-site vacancy (Pt-^Pt(V)) and Pt 
to Co-site vacancy (Pt— >Co(V)). In Fig. 1 we represent 
the two possible jumps of a Pt atom toward a vacancy 
situated either in the same plane (Pt— ►Pt(V)) or in a 
neighbouring plane ( Pt— >Co(V)). In the case of a relaxed 
lattice, we have used a Tight Binding-Quenched Molecu- 
lar Dynamics-Second Moment Approximation (TB-QMD- 
SMA), which allows us to determine the equilibrium struc- 
ture of a system with a finite number of particles at T = 
K, by integrating the equation of motion |37|. The quench- 
ing procedure, in which the velocity vi of an atom i is 
cancelled when the product Fi(t)vi(t) is negative, leads 
to the minimisation of the potential energy at if |38| . 
Fi being the force acting on the atom i, calculated in the 
extended tight-binding formalism from the total energy 
|34| . The positions were integrated by means of the Verlet 
algorithm [21] 

For all atomic jumps, the total energy follows a smooth 
cosine curve, with a maximum close to half-way between 
the initial and destination lattice sites. In the case of the 
cross jumps, Co^Pt(V) and Pt— >Co(V), a shift in the to- 
tal energy, associated to the formation of an anti-site de- 
fect, was observed at the final position. The saddle-point 
energies have been deduced from the evolution of the total 
energy for the different kinds of atomic jumps. The results 
obtained with and without relaxation are reported in Ta- 
ble III. The saddle-point energies are denoted e R and e for 
the relaxed and the non-relaxed lattice respectively. The 
amplitude of the saddle-point energies may be explained 
by the size effect and the atoms forming the barrier of first- 
nearest neighbours that the jumping atom has to overcome 
for completing the jump. For instance, the largest value is 
obtained in the case of Pt— >Pt(V) jump, when the bigger 
Pt atom is involved in the jump. In the case of the cross 
jumps, Co— ►Pt(V) and Pt— >Co(V), the atomic structure 
of the barrier to overcome being the same for the two 
jumps, the corresponding saddle-point energies are very 
close. Furthermore, the lattice being softer when relaxed, 
the saddle-point energies, are as expected, smaller when 
the relaxation is taken into account. 



3 Monte Carlo simulations 
3.1 Simulation method 

MC simulations have been established as a useful tool for 
studying order-order and order-disorder relaxation kinet- 



ics in intermetallics |4()H41U42j . In contrast to earlier stud- 
ies, based on effective pair interaction energies expressed 
within a simple Ising-Hamiltonian, we go in the present 
work beyond the pair approximation by implementing the 
many-body potential deduced from the TB-SMA and tak- 
ing into account the saddle-point energies calculated in the 
previous section to simulate the disordering process of a 
perfectly ordered Llo-CoPt structure. The purpose of the 
simulation is to prove whether the interaction model is 
appropriate and to get an estimation of the vacancy mi- 
gration energy which is an important energetic parameter 
that drives diffusion and ordering process in intermetallics. 

In view of recent neutron diffuse scattering measure- 
ment in equiatomic CoPt |35| which show a symmetric dis- 
tribution of the diffuse intensity around the 100 and equiv- 
alent points in the reciprocal lattice, signature of a highly 
stable Llo phase with very small static displacements, we 
have chosen to carry on the MC simulations using a rigid 
lattice. We have used a model based on the vacancy-atom 
jump mechanism between nearest-neighbour sites, which 
is the realistic microscopic process in dense phases [2*9] . 
The simulation box contains 32 3 Llo fcc-based cells with 
linear periodic boundary conditions. The simulation starts 
with a perfect Llo ordered crystal in which one of the 
two sublattices (sublattice a) is occupied by Co atoms 
and the other (sublattice (3) by Pt atoms. To not affect 
the static properties and avoid interaction effects, a sin- 
gle vacancy is introduced at random in the crystal. The 
elementary MC step is the following: one of the vacancy 
neighbours (Co or Pt atom) is randomly chosen, the en- 
ergy balance AE of the atom- vacancy exchange before and 
after the jump is evaluated. The jump is performed if the 
Glauber probability P{AE) = P G {AE)exp{-e/k B T) is 
larger than a random number between and 1. P G ( AE) = 
[1 + exp(AE I fceT)] -1 is the Glauber probability 03] and 
e is the saddle-point energy for the corresponding jump. 
This corresponds to averaging the result over a large num- 
ber of reversal jump attempts, the sum of the probabilities 
of the jump and its reversal being equal to 1. 

The configuration of the system was analysed, at reg- 
ular time intervals, by calculating the LRO parameter 
rj = 2{2Ng a - N C o)/{N sitcs - 1), where Ng a is the num- 
ber of Co atoms on the a sublattice and Nco the total 
number of Co atoms. The time scale being the number 
of jump attempts. For each temperature, the evolution 
of rj is followed until the system reaches equilibrium. We 
have chosen to stop when the simulation time was at least 
longer than 5 times the relaxation time of the system. 



3.2 Results and discussion 

Isothermal relaxations of 77 have been recorded in the tem- 
perature range from 600 to 980 K. An example of kinetics 
relaxation is shown in Fig. 2. According to the path prob- 
ability method |44j . corroborated by many experimental 
results 1"> l(i 1711 1811 the isothermal relaxation of the 
LRO parameter in intermetallics is well fitted with two 
exponentials, yielding a long and a short relaxation times, 
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corresponding to the slow and the fast processes, respec- 
tively. This result is well established in LI2 phase for which 
a detailed study has shown that the fast process is related 
to the formation of the nearest-neighbour antisite pairs, 
whereas the slow one is related to the uncoupling of these 
antisite pairs Due to the difference in structure, the 
process must be different in the Llo phase and the inter- 
pretation is still under investigation. However, recent MC 
simulations of the LRO relaxation in Llo-FePd showed 
that the fast process is highly predominant (up to 90%) 
below the order-disorder transition 50 . In view of these 
results and considering the temperature range of our sim- 
ulations, we have chosen to fit the kinetics using a sin- 
gle exponential yielding a single relaxation time r which 
fulfils an Arrhcnius law with a positive migration energy 
Em = 0.73 ± 0.15 eV (Fig. 3). The relaxation times are 
typically in the order of 10 12 MC steps, significantly higher 
than the relaxation times obtained without saddle-point 
energies. The data points on the left side of Fig. 3 , plotted 
in open circles, show a clear departure from the straight 
line, they are synonym of a critical slowing down close to 
the order-disorder transition and have been excluded from 
the linear regression. 

The vacancy formation energy has been measured 
in pure Co and pure Pt, it has been found equal to 1.38 eV 
and 1.2 eV for Co and Pt respectively [3]. Assuming that 
the activation energy E^ is the sum of Em and Ep, the 
calculated value of E^ obtained considering the simulated 
value of Em and an average value of Ei? compares well to 
the activation energy of 2.52 eV, measured in equiatomic 
CoPt by means of resistivity measurements |28| . Normal 
modes of vibration in equiatomic CoPt have been recently 
measured by inelastic neutron scattering |52|. The phonon 
density of states has been used to calculate the migra- 
tion energy using Schober's model. A value of 0.85 eV 
was found in the fec-disordered state at 1120 K. Despite 
of the difference in the state of order between the simu- 
lated CoPt system and the measured one, the qualitative 
agreement is quite satisfactory. In our knowledge no fur- 
ther measurement of the migration energy in equiatomic 
CoPt are available. Nevertheless, qualitative comparison 
to CoPt3, in which the migration energy was calculated 
using both Schober's and Flynn's models, can be made. 

Schober's model initially developed for bec and fee 
pure metals ^H] has been extended to the A3B compounds 
with Ll 2 structure JT] and applied to CoPt3 - Aver- 
aging the values of the migration energy over the differ- 
ent kinds of atomic jumps, we get a migration energy of 
1.3 eV at 300 K. Flynn's model was applied to CoPt 3 
using the elastic constants deduced from the slope of the 
phonon dispersion curves at the center of the Brillouin 
zone jSSj- A migration energy of 0.95 eV was found in 
the fec-disordered state at 1060 K. The difference in the 
amplitude between the calculated value of the migration 
energy in CoPt and the values deduced from Schober's 
and Flynn's models in CoPt3 might be explained by the 
atomic mass and size effects. As expected, an increase in 
the migration energy with the atomic density is observed 



between the Llo and LI2 phases of Co-Pt system. The 
atomic density is nearly 25% lower in CoPt than in CoPt3. 

The variation of the equilibrium LRO parameter r\ eq as 
a function of temperature is shown in Fig. 4. The order- 
disorder transition takes place between 930 and 935 K. 
The transition region has been crossed with 5 K tem- 
perature steps. It should be noted that the value of this 
temperature is close to the value obtained in the Llo com- 
pounds using phcnomcnological pair interaction energies 
[5~1"| but it is 15% lower than the experimental value of 
the order-disorder transition temperature T c in equiatomic 
CoPt (1110 K) Nevetheless, the general agreement 
between the simulated and the experimental T c is quite 
satisfactory. In fact, the MC model, which assumes tem- 
perature independent potentials and does not take into 
account anti-phase domains, only gives an estimate of T c . 
Indeed, recent x-ray diffraction and transmission electron 
microscopy measurements in CoPt showed that the order- 
ing transformation involves formation of antiphase bound- 
aries and twin bands |17|. 

4 Conclusion 

An approach for determining the parameters of a many- 
body potential in equiatomic CoPt using the TB-SMA 
has been presented. The potential was used to determine 
the saddle-point energies for the different kinds of nearest- 
neighbour atom- vacancy jumps in the case of a rigid and a 
relaxed lattice. The calculated energy barriers have been 
used together with the many-body potential to simulate 
the LRO relaxation in Llo-CoPt structure using a MC 
technique. A vacancy migration energy of 0.73 ± 0.15 eV 
has been deduced from the Arrhcnius plot of the relax- 
ation times. A lack in experimental data of the migration 
energy in the ordered state of the equiatomic CoPt made 
possible only qualitative comparisons to either the disor- 
dered state of CoPt data, when available, or to CoPt 3 
data. Finally, the order-disorder transition temperature 
has been determined from the isothermal relaxation of the 
LRO parameter. The simulated value is in a satisfactory 
agreement with experiment. 
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Fig. 1. The two possible atomic jumps, in-plane and out-' 
plane, of the Pt atom. 
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Fig. 2. Isothermal relaxation of r\ obtained for T = 880 K 
(black circles) and its simulation using a single exponential 
(white line). 
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Fig. 3. Arrhenius plot of the relaxation times. The linear re- 
gression gives the migration energy Em = 0.73 ±0.15 eV. The 
open circles show the data points which have been excluded 
from the linear regression. 
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Fig. 4. Temperature variation of the equilibrium LRO param- 
eter r/ eq . 
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Table 1. Potential parameters of inter-atomic interactions de- 
duced from the TB-SMA. The interaction parameters between 
Co-Co and Pt-Pt are taken from the literature |36) 



Interaction 


A (eV) 




P 


q 


Pt-Pt 


0.242 


2.506 


11.14 


3.68 


Co- Co 


0.189 


1.907 


8.8 


2.96 


Pt-Co 


0.175 


2.115 


9.412 


2.812 



Table 2. Reduced lattice parameters, calculated (E ca i) and 
experimental (E exp ) formation energies of the ordered Llo and 
disordered Al phases. In Llo phase, the results obtained by 
varying a, c, and both a and c simultaneously are shown. We 
have used ao = 3.806 A, Co = 3.684 A for the Llo phase and 
a = 3.751 A for the Al phase 





E ca i (eV) 


E exp (cV) 


a/ao 


c/c 


Llo (varying a) 


-0.133 


-0.14 


0.99 




Llo (varying c) 


-0.142 


-0.14 




0.97 


Llo (varying a and c) 


-0.138 


-0.14 


0.995 


0.995 


Ai (varying a) 


-0.086 


-0.10 


1 





Table 3. Calculated saddle-point energies in eV for the differ- 
ent kinds A— >B(V) jumps of an atom A on a vacancy site B in 
the case of a rigid (e) and a relaxed lattice (e R ). 





Co^Co(V) 


Co^Pt(V) 


Pt-»Pt(V) 


Pt^Co(V) 


e 


0.49 


0.41 


0.54 


0.39 


e R 


0.39 


0.34 


0.42 


0.31 



